Probability-based comparison of quantum states 



O 
(N 

in 



Sergey N. Filippov^ and Mario Ziman^'^ 
^Moscow Institute of Physics and Technology, Moscow Region, Russia 
'^Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland 
'^Institute of Physics, Slovak Academy of Sciences, Bratislava, Slovakia 
(Dated: February 5, 2012) 

We address the following state comparison problem; is it possible to design an experiment enabling 
us to unambiguously decide (based on the observed outcome statistics) on the sameness or difference 
of two unknown state preparations without revealing a complete information about the states? We 
find that the claim "the same" can never be concluded without any doubts unless the information is 
complete. Moreover, we prove that a universal comparison (that perfectly distinguishes all states) 
also requires the complete information about the states. Nevertheless, for some measurements, the 
probability distribution of outcomes still allows one to make an unambiguous conclusion on the 
difference between the states even in the case of incomplete information. We analyze an efficiency 
of such a comparison of qudit states when it is based on the SWAP-measurement. For qubit states, 
we consider in detail the performance of special families of two-valued measurements enabling us 
to successfully compare at most half of the pairs of states. We conclude with a surprisingly simple 
example of an almost universal comparison measurement in any dimension. 

PACS numbers: 03.67.-a, 03.65.Wj 
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INTRODUCTION 



The exponential scaling of the number of parameters 
describing multipartite quantum systems stands behind 
the potential power of quantum information processing. 
However, the same feature makes a complete character- 
ization (tomography) of unknown quantum devices in- 
tractable. Therefore, it is of practical interest to under- 
stand which properties of physical systems require the 
full tomography for their determination and for which of 
them such a complete knowledge is redundant. In this 
paper we analyze the resources needed for a comparison 
of quantum states. Suppose given a pair of quantum 
systems in unknown states. The question is what ex- 
periments (if any) are capable to reveal with certainty 
either the difference between the states, or confirm their 
sameness as long as the the probability distribution of 
measurement outcomes is identified. 

By the very nature of quantum theory the events we 
observe in quantum experiments are random. That is, 
both quantum predictions and quantum conclusions are 
naturally formulated in terms of probabilities and uncer- 
tainty. Therefore it is surprising that there are (very spe- 
cific) situations (including special instances of the com- 
parison problem) in which individual clicks enable us to 
make a nontrivial unambiguous prediction, or conclusion. 
For example, if we are given a promise that the states are 
pure, then (with a nonzero probability) the difference of 
states can be confirmed unambiguously from a single ex- 
perimental click [T, '51 . This result can be also generalized 
to the comparison of many pure states , the compar- 
ison of ensembles of pure states @ , and the comparison 
of some pure continuous-variable states 0, Q (see also 
the review [9,])- Unfortunately, such single-shot (non- 
statistical ) co mparison strategy fails for general mixed 
states [1, [lo| . The reason is simple. The probability 



of any outcome is strictly nonvanishing providing that 
a bipartite system is in the completely mixed state, for 
which the subsystems are in the same state. That is, for 
any outcome there is a situation in which the systems 
are the same, hence the difference cannot be concluded 
unambiguously. In such case any error-free conclusions 
need to be based on the observed probabilities of out- 
comes. Probability-based comparison strategies were not 
considered previously, so we fill this gap in the present 
paper. 

Trivially, if the experimentally measured probabilities 
provide a complete information on quantum states of 
both systems individually, then they contain also all the 
information needed for the comparison. The question of 
our interest is whether the complete tomography is neces- 
sary. Our main goal is to design a comparison experiment 
providing as little redundant information as possible. 

In Sec. [ni we introduce the necessary mathematical 
notation and formulate the problem. In Sec. IIIIl we ad- 
dress the existence of a universal comparison measure- 
ment. Sec. IIVI investigates the comparison performance 
of two-outcome measurements. An example of an almost 
universal three-outcome comparison measurement is pre- 
sented in Sec. IVland conclusion is the content of Sec. lVIl 



II. 



PROBLEM FORMULATION 



Any quantum state is associated with the density oper- 
ator Q G S{T-L) such that g> O and tr[g] — 1. Hereafter, 
S{H) stands for the set of all states of a system associ- 
ated with the Hilbert space H. The statistical features 
of quantum measurements are fully captured by means 
of a positive operator-valued measure (POVM) that is 
a collection E of positive operators (acting on 7i and 
called effects) Ei,...,En summing up to the identity. 
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FIG. 1: Illustration of probability-based comparison. If the 
observed probability distribution belongs to \ , then 
states g and ^ are for sure different. 

i.e. J2]j=i ^3 — ^- For each state g € 5(H) the measure- 
ment E assigns a probability distribution = Pe, 
where pj — tr[Ejg] > and X]j=iPi — 1- 

Let us now move on to the set of bipartite factorized 
states 5fac = {e «> f : £>, C S Sin)} C S{n (g) n), where 
the parties g and ^ are the states to be compared. For a 
fixed measurement E we can ask how much information 
it reveals concerning the comparison of the subsystems. 

Denote by S'^ the subset of twin-identical states, i.e. 
S+ = {rj^Tj : rj £ S{n)} C S{n (g) n). Similarly, let 
us denote by the subset of non-identical states, i.e. 
S~ = {g®^ : g,^{^ g) £ S{n)} C S{H®H). Obviously, 
iSfac = S^yjS~ . The goal of comparison is then to distin- 
guish between sets of states S'^ and S" . This goal can be 
achieved in our approach by considering two sets of prob- 
ability distributions = {p:pj = tr[£jw], cj g 5*}. In 
other words, since the measurement E performs the map- 
ping one can unambiguously conclude that a 
bipartite state w belongs to the set if the observed 
probability distribution p^ £ \ (see Fig. [1]). 

For a fixed POVM E on H ig) "H we may introduce the 
following quantities: 

DE{g(^^,S+) = inf V|tr[^;,(^.®^-?7®r,)]|,(l) 

3 

DE{'n®'>hS^) = inf \iT[EAg®^- ri®'q)]\.{2) 

While D^{g®£^,S^) quantifies how different the states 
g and ^ are (with respect to measurement E), the value 
of 0^(7] (E) f]j'S~) tells us to which extent the equivalence 
of twin-identical states can be confirmed. 

Before we proceed further let us make one important 
observation: for all e > and any state 77 g) 77 there exists 
a state g^ ^ such that |tr[i?(77 (E) rj — g (E) < e for any 
POVM effect E. In other words, in order to conclude 
that the states are the same no uncertainty in the spec- 
ification of the probabilities pe{'^) — tr[i?w] is allowed. 
Such an infinite precision is practically not achievable, 
however, for our purposes we will assume the probabil- 
ities are specified exactly. The proof of the statement 
above is relatively straightforward. Let us set g = -q and 
i = (l-D'y+^-f, i-e. ri<^ri-g<^^ = ^■q(E{'q-\l). Since 



|tr[£;X]| < maxseEtr[|£:X|] < tr[|X|], it follows that 

\tr[E{7^®ii- g®m<^\^'^[\v®{n--J)\]<^- (3) 

In the last inequality we used the fact that the trace 
distance of states is bounded from above by one. 

Our observation implies that for any measurement E 
we have !?£(?? ® Vi^ ) = 0- This seems to be in con- 
tradiction with measurements which provide us with a 
complete information on the states of individual systems. 
We will refer to such measurements as locally informa- 
tionally complete (LIC) measurements. Clearly, in case 
of an ideal LIC measurement the sameness can be ver- 
ified. Where is the problem? Topologically, in the set 
of factorized states 5fac with the trace-distance metrics, 
the subset S'^ is closed and does not contain any interior 
point (therefore the distance ([2]) vanishes), however, it 
does not mean that the subset 5+ is empty. 

The vanishing value of the considered distance is not 
completely relevant if one thinks about ideal error-free 
experiments. In practise, the experimental noise is un- 
avoidable, hence, from the practical point of view a con- 
clusion on the sameness of states can never be error-free. 



III. UNIVERSAL COMPARISON 
MEASUREMENT 

We say the measurement E implements the comparison 
whenever DE{gE)£,,S^) > for some pairs g,£,. The 
state comparison measurement E is universal if D£{g (E 
^, 5*+) > for all g, ^(^ g). This is, for instance, achieved 
in case of ideal LIC measurements even though the value 
of D^ig E) ^,S~^) can be arbitrarily small, but always 
strictly positive. As before, this situation is in practice 
not very realistic, because any error in identification of 
the outcome probabilities makes the conclusions (in some 
cases of g and ^) ambiguous. However, assuming the 
infinite precision in specification of the probabilities, the 
universality can be achieved and in what follows we will 
assume that probabilities are identified perfectly. The 
potential errors can be viewed as modifications of the sets 
and S~ we are aiming to distinguish. Nevertheless, 
our goal is to analyze the ideal case. 

Let us now demonstrate that a universal comparison 
can be implemented if and only if the measurement is 
LIC. 

To start with, we are reminded that any POVM E with 
effects Ej linearly maps a state oj £ S{H E) H) into the 
probability vector p = (pi,p2, • • where pj = tT[EjUj]. 
For LIC measurements the induced mapping p g) ^ 1— ?■ 7? 
is bijective. That proves the sufficiency. To prove the 
necessity let us assume the converse, i.e. suppose the 
measurement E is not an LIC measurement but imple- 
ments a universal comparison. Since E is not LIC, the 
probability assignment g (g ^ H> pe is injective. Let us 
denote by 11^ and 11 the images of and iSfac under 
some LIC measurement, respectively, and by Pe denote 
the image of iSfac under measurement E. Clearly, the 
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relation between Tr{g(S) £,) and 'Pe{q'^ £,) is linear and in- 
jective, i.e. there exist probability vectors tti S 11 and 
TT 2 e n transformed into the same probability distribu- 
tion pe{q ^ € ^E- The distributions ifj transformed 
into the same probability vector pe{q ® £,) span a linear 
subspace (hyperplane) Hg(g,^ in the linear span of 11. 

Consider now a non-pure state rj ® -q, then its image 
'rf{rj ® rj) is an interior point of 11 on the probability sim- 
plex. There exists Co > such that for all < e < eo the 
neighborhood Oe(7r (77(g) 77)) belongs to 11 (on the simplex). 
Moreover, the intersection Oe(7f(?7 (g) 77)) n cannot 
be a subset of 11+ only, because 11+ does not contain 
any interior point on the simplex (if it did, the distance 
£'e(?7 €5 Vt'S ) would not vanish for all states 77). Thus, 
Oe(7r(77(g)7/))niJ^0,,nn~ is not empty and contains points 
of the form n{g(^^) such that g ^ As both 'k{t]®t]) and 
■K{g®C) belong to i?,,®,,, we have pe(^ <8) |) = Pe(?7 <8) 77), 
i.e. E is not a universal comparison measurement. This 
contradiction concludes the proof of the necessity. 

Let us summarize two main conclusions: 

(i) In any locally informationally incomplete measure- 
ment the sameness of states cannot be confirmed. 

(ii) Universal comparison (concluding universally and 
unambiguously the difference of states) requires a locally 
informationally complete measurement. 

A question that remains open is how to evaluate the 
overall performance of (universal, or non-universal) com- 
parison experiments. There are several options. We can 
use the volume of the subset S~^^p of states in S~ that 
can be successfully compared, or the average value of 
DE{g®£,,S^) with respect to some measure on the state 
space. In particular, these quantities read 

l^complE = jj^ fi{dg)fi(dOh{DE{g<E>^,S+)), (4) 

(De) - JJ^ ^,{dg)^lidODE{g(E>^,S+), (5) 

where h{x) is the Heaviside function and iJ-{dg) — fi{d^) 
is a measure on the state space of individual subsystems. 
Quite common choices for the measure /i on density oper- 
ators are the ones induced by metrics, namely, by Bures 
distance and Hilbert-Schmidt distance (see, e.g., fl^lTsj 
and references therein). Let us stress that l^^^^mplE = 1 
does not imply the comparison is universal, because there 
can be a set of measure zero for which DE{g^(,, S'+) = 0. 
In such case we say that the comparison measurement 
is almost universal. It is of great interest to investigate 
whether there exist some almost universal comparison ex- 
periments and, in particular, how many outcomes such 
measurements require. 



IV. TWO- VALUED COMPARISON 
EXPERIMENTS 

Let us start our investigation with the simplest case of 
two- valued POVMs described by the effects E,I — E. In 



such case 

DE{g^^,S+) = 2DE{g®^,S+), (6) 

where Dsig ^,S^) = min^^^g5+ \tTE{g ® ^ - 77 (g) 77)]. 
Such measurements cannot be LIC, thus, they are nec- 
essarily non-universal comparators. Nevertheless it is of 
practical interest to understand how good is their com- 
parison performance. 

A. SWAP-based comparison 

As we have already mentioned in Sec. HI if we restrict 
ourselves only to pure states, then there exists a strategy 
to perform an unambiguous comparison (via the SWAP 
measurement). In such an approach, the sameness of the 
states cannot be concluded and this is related to the ab- 
sence of the universal NOT-operation fTl|. However, the 
strategy (if successful) can reveal the difference between 
the states in a single shot, hence, no collection of statis- 
tics is needed. 

The key observation for such a conventional strategy 
is that the support of twin-identical pure states spans 
only the symmetric subspace oi % ® %. Suppose pro- 
jections -Esym , -E'asym onto Symmetric and antisymmetric 
subspaces of "H g) "H. Since i?sym + E^sym = I they form 
a two- valued POVM Eswap- Let us note that i?sym = 
\{I + S), E^sym = i(/-5'), where S is the SWAP opera- 
tor acting as S{\ip^(p)) — \(p(^tp) for all \tp), e H. It 
is straightforward to see that for any twin-identical pure 
state \(p(d(p) one has tr[£'asym|</5(g<i2)((^g)<p|] = 0, however 

tr[Sasym|^'»^)(<^®^|] = i(l-|(^|^)H > if \^) ^ \if) . 

Therefore, recording an outcome -Easym allows us to un- 
ambiguously conclude that the states are different. No 
statistics is needed. 

Let us see how this strategy works in the case of general 
mixed states. A direct calculation yields 

ftym = tliE,yn,g g) = ^(1 + tr[f?^]) , (7) 
Pasym = tr [i;asym g g) ^] = 5 ( 1 - tr [g^] ) , (8) 

where we used the identity tr[S'g(g^] = tr[p^]. The purity 
of states tr[7;^] is bounded from below by 1/d, where d — 
dimH. Thus, 

Psym(e0C) e [il) = Ps;m, (9) 
ftym(?7®^)e (10) 

where is the image of 5+ under the POVM ef- 

fect Esym- It follows that by measuring the probability 
Psym < {d+l)/2d we can with certainty conclude that the 
states are different. In particular, DgwApig ® fi'5+) = 
max{0, i - tr[g(\}. 

Suppose g — i(/-|-r-A) and = i(/-|-k-A), where A = 
(Ai, . . . , A^2_i) is a vector formed of traceless hermitian 
operators Aj such that tr[AjAfe] = dSjk, and r, k G M'' ~^ 
are Bloch-like vectors. Using this notation, we find that 
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FIG. 2: Body _B(xo), i-e. the region of parameters (xi, >tz) when (|12|) is a true POVM effect. Parameter j<o takes values i, 
|, |, and I for figures from left to right. The union Uj<(jg[o,i]-B(>fo) is a rhombododecahedron and is depicted by solid lines. 
POVM effects -Easym, Saym, i?z± are vertices and POVM effects -Eiy± are face centers of this convex polytope. 



-DswAp(g ® '5^) > only if r • k < 0. That is, for each 
fixed r(^ 0) the set of successfully comparable Bloch-like 
vectors k is given by an intersection of the hemisphere 
with the state space. 

Let us stress that for qubits [d ~ 2) it is exactly the 
hemisphere, thus, l^^^mplswAP = \ (see the next para- 
graph). In other words, a difference of states from the 
same hemisphere is not detected in the SWAP measure- 
ment. This implies that the approximate universality is 
lost. 

Due to unitary invariance of the measure \x we can 
always treat one of the states in De{q ® ^,5+) as di- 
agonal, say Q. Then tr[gi^] — Yl'j=i Qjj^jj ^^.d the in- 
tegration area of |iS^jj^p|swAP is split into d\ subsets la- 
belled by the permutation of the labels ji, . . . , jd iden- 
tifying the ordering gj^j-^ > ■ ■ ■ > Qj^jd of eigenvalues 
of g. The (normalized) volume of each of these sub- 
sets is ^. If and ^ are from mutually opposite sub- 
sets (labelled as ji, . . . ,jd and jd, ■ ■ ■ respectively), 
then tr[g(] ~ Qjj^jj — 2 meaning that such pairs 

g and can be successfully compared. Therefore, we 
found a lower bound |5(^j„p|swAP > ^- If C are from 
the same subset, then tr[g£] — X]j=i Qjj^jj — 2' hence 
the contribution to |5(^„jp|swAP is vanishing and we can 
bound the fraction of comparable states also from above, 
I'^complswAP < 1 — 3T- For qubits both bounds coincide 
and give the value |5(;omplswAP = ^■ 

B. Qubit "diagonal" comparison experiments 

In the previous section we have shown that the SWAP- 
based comparison enables us (in case of qubits) to suc- 
cessfully detect the difference for half (up to a set of zero 
measure) of the pairs g(E)^- Can one do better with some 
other two- valued measurement? 

Let us stress that, for a given effect E, the states g 
comparable with a fixed state ^ form a cut of the Bloch 
sphere. Therefore, if |5(5^^p|swAP > ^, then the com- 



plete mixture must be comparable with a subset of a non- 
vanishing measure. Note that for the SWAP-based com- 
parison complete mixture cannot be conclusively com- 
pared with any other state. 

A general two-qubit effect takes the form 

where we use the notation ctq = I and cti , (72 , denote 
the Pauli operators (Ta; , tij, , , respectively. Real coeffi- 
cients Elm read — jtT[Eai (E) Urn\ and satisfy the the 
constraints O < E < I . 

For the sake of simplicity let us consider only the diag- 
onal case, i.e. we will assume sim = SimXm, with >Cm € K 
for all m = 0, . . . , 3. Then 

/ >fO + >f3 Xi ~ X2 \ 

_ >Cq - >C3 + I 

'^'''^ ~ I >fi + >f2 ><0 - 

\ Xl — X2 >io + X3 / 

(12) 

Applying the positivity constraints on E of this form we 
obtain the following conditions; 

^ — — X2 — ^ 1, (13) 

< >co — xi + + ^3 l£ 1, (14) 

< >fo + >fi ^ >'^2 + < 1, (15) 

< >iro + >fi + >f2 ^ ^3 !i 1- (16) 

System of inequalities p^ -([T5 | has a nontrivial solu- 
tion whenever < xq < 1. Indeed, if xq is fixed, then 
each two-sided inequality of the system determines a ge- 
ometric fiber between two planes in the reference frame 
(>fi, X2, with the distance between the planes being 
equal to It < >cq < j, then four fibers intersect 
to yield a tetrahedron. The intersection becomes a trun- 
cated tetrahedron if i < >fo < | and finally transforms 
into an octahedron for the case xq — ^. li xq > ^ 
then the solution is a body obtained by the inversion 
(>fi, >f2, X3) — {—xi, —X2, —X3) of the body labelled by 
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parameter (1 — xtq). For instance, if | < >fo < 1: then 
the intersection is a tetrahedron inverted with respect to 
that in case < < j. Given hq we will refer to the 
intersection as body B{>to) (see Fig. 

The probabilities of the measurement outcome, corre- 
sponding to the POVM effect i?diag, can be readily calcu- 
lated for the non-identical (different) and twin-identical 
(same) states and read 



Pdiff 

Psaiac 



Pdiag{g <^£.) = ^0 + J2rn=l ^^mJ'mSm , (17) 
-3 ,2 



(18) 



where we used g ~ i(/ -I- r • cr), ^ = i(/ -|- s • cr) and 
•q = i(/-|- h • cr). Using the normalization constraints for 
Bloch vectors r, s, and h, we obtain from Eqs. PT )) - P^ 
that probabilities pdm and Psamc satisfy the relations 



Pdiff e [XO — >^max, ><0 + ><maxj 



Ps. 



(19) 
(20) 



respectively, where we used >!:inax — niax{ \xi\, \xi\^ | >?3 1 } , 
X- = min{0, >fi, >f2, >?3}, and x+ = max{0, >(ri, >f2, >^3}- 
Clearly, J^^iag ^dlag' hence, there is a two- valued 
POVM that allows making a nontrivial conclusion on the 
difference of some states. Notice that the case x-q = |, 
xi = X2 — x-i — J gives the SWAP-based comparison 
measurement, for which whenever the measured proba- 
bility Pdiag satisfies Pdiag < I J then the states g and ^ are 
unambiguously different. 



1. Fixed pair comparison 

Surprisingly, there are pairs of states g and ^ such that 
no measurement E of the form (fT2| can reveal their dif- 
ference. A direct calculation gives that the difference for 
a pair of states g and ^ can be concluded in the diagonal 
comparison experiment if 



-Ddiag(£i<8'€,'5+) = min 

|h|<l 



X/m=l Xmi^mkm ^m) 



> 0. 



(21) 

Suppose that rmkm is nonnegative for all to. Setting 
hm = V^^mkm the distance (|2ip is vanishing for arbitrary 



measurement of the considered diagonal form. Let us 
stress that the requirement of positivity of rmkm for all 
TO means that signs of the Bloch vector components co- 
incide, hence, k and r belong to the same octant of the 
Bloch ball. Let us stress, however, that the octants de- 
pend on the choice of the axes (Pauli operators), and for 
a given pair of states we can always fix the coordinate 
system in such a way that they belong to two different 
octants. The only exceptions are collinear vectors k and 
r = ck for c > 0. In fact, a pair of parallel Bloch vectors 
(pointing into the same direction) is indistinguishable by 
any diagonal measurement irrelevant of the choice of the 
coordinate system. In particular, it follows that none 
of these measurements is capable to distinguish (in the 



comparison sense) the complete mixture g = from any 

other state, because Pdis{^I ^0 = Psamd^I ^ ^I) — xq. 

It is natural to ask for which pairs g, ^ their difference 
can be identified by a suitably selected E of the consid- 
ered diagonal form and whether there are some "non- 
diagonal" measurements enabling us to compare a pair 
of states containing the complete mixture. 

In order to get an insight into the power of diagonal 
measurements, let us assume that k„i > 0, m = 1,2,3 
and fix k. Define a new vector K = {>ciki, >C2k2, x^k^). 
If K • r > 0, then we can find h such that K • r = 
Lm^i^m^m, hence Ddiag = 0. If K • r < 0, then 
i?diag = |K • r| + J2L=i ^™^m > for any h ^ 0. There- 
fore, the minimum is achieved for rj ® rj = i/ (8) ^/ and 
the distance reads 

^--(^^^"5") = {|K.r| I'^rwis-e."' ^^2) 

In other words, the considered diagonal measurement 
E enables us to verify the difference of g and ^ for all g 
satisfying the inequality K • r < 0. The condition K • r = 
determines a plane containing the complete mixture 
(center of the Bloch ball), hence, for any measurement of 
the considered type and any state ^ the set of successfully 
comparable states g is exactly a hemisphere of the Bloch 
ball. Fig. [3] illustrates this situation for the following 
choices of the diagonal measurements (POVMs): 



(g)a„)}; (23) 

Cm ® dm), 

(24) 



EsWAP 


- lE 


= K3 


I ® 


^ + Em=l 




-^asym 


= \{I^^ 


5/- 


2^m=l 




= {Exy^ 


= i(3 


I® 


I + YL^i 




E 




I- 


X]m=l ^rn 






= i(/« 


I± 


0-3 (g) era)}. 



(25) 



In particular, for EswAP the "comparable hemisphere" 
is orthogonal to the vector k. For E^-y the "comparable" 
hemisphere is orthogonal to the vector ky — (fci,fc2,0) 
being a projection of k onto the xy plane. Finally, for E^ 
any state from the northern hemisphere is "comparable" 
with any state from the southern hemisphere. 

It is worth noting that we have restricted ourselves to 
a specific form of POVM effects (|12p. However, even 
for such a simplified problem the solution looks rather 
sophisticated. 



2. Average performance 

The fact that for any given diagonal two-valued mea- 
surement the states within the same octant are not com- 
parable means that none of them is universal neither in 
an approximative way. Nevertheless, it is of interest to 
understand which of them perform better than the oth- 
ers and which do not perform at all. In particular, we 
are interested in the answer to the following question: 
How many qubit states g, ^ can be distinguished? As it 
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(a) (b) (c) 




FIG. 3: States g in the Bloch ball (determined by vectors r) which can be distinguished from a fixed state ^ (given by vector 
k) by using measurement Eswap (a), Exy (b), and Ez (c). 



is briefly outlined in Sec. Illlj to answer this question it 
is necessary to introduce some measure /z on the state 
space. Once it is introduced, we can evaluate the quan- 
tities IS", 



write 



comp I 



(relative volume of the successfully com- 
parable states, i.e. the comparison universality factor) 
and (-De) (average distance, i.e. the comparison quality 
factor) . 

In contrast with pure states, for density operators there 
exist many equivalently well-motivated measures (see, 
e.g., and references therein). We will employ two 

most commonly used ones, namely, the Hilbert-Schmidt 
measure /ljhs and the Bures measure 



msidg) = 

fiBidg) = 



47r 



r sin 9 dr dd dip , 



: dr dO dip , 



(26) 
(27) 



where we used the following parametrization of Bloch 
vectors: r — {r cos ip sin , r sin ip sin 9 , r cos 0) with r G 
[0, 1], £ [0, tt], and p £ [0, 2tt]. Both these measures are 
spherically symmetric and the former one corresponds to 
the uniform coverage of the entire Bloch ball [l2j, i.e. 
the Hilbert-Schmidt measure /ijjs (T) of any compact set 
T C 5(7^2) equals the geometrical volume Jg(^^-^^rp d'^r of 
the corresponding body inside the Bloch ball divided by 
47r/3. The Bures measure ((?7l) ascribes higher weights 
to the states with higher purity (that are closer to the 
surface of the Bloch ball). 

For the calculation purposes it is convenient to intro- 
duce the following (relative) density of states: 



lim — — 

Ap^O Ap 



ti[Eg<»^]<£[p-p+Ap] 



fi{dg)fi{d^) , 



(28) 

whose physical meaning is that Nf^^{p)l\p equals the 
fraction of pairs g® ^ resulting in the measurement out- 
come probability within the region [p,p -f Ap] for the 
effect E. Using the introduced function, we can readily 



1^. 
(Ce) = 2 



_ diff 

comp I E / ■'■'^ E,fi 



N%^l{p)dp, (29) 
\p-Po\N'il{p)dp, (30) 



where stands for the image of under the action 
of POVM effect E and pq is simultaneously the frontier 
point of and the inner point of (if there are two 
such points, then (De) is a sum of two integrals ([50)1 in 
the corresponding regions of variable p) . 

In what follows we will analyze the three examples 
from the previous section (comparison measurements 
Eswap, ^xy, and E^) and compare their performance. 
The associated densities are depicted in Fig. 0] and explic- 
itly written in the Appendix. We focus on these POVMs 
because they represent three different types of boundary 
(extremal in case of EgwAP and E^) points of the set of 
diagonal measurements (see Fig. ^ . Other diagonal mea- 
surements will exhibit intermediate behavior with respect 
to these three. 

We have numerically tested that for diagonal mea- 
surements the relative volume of the set of compara- 
ble states |5(^jj^p|e < 5 irrelevant of the measure used. 
The considered three examples saturate this value, i.e. 
l'^complEswAP,ccy,2 — \- In fact, there exist many measure- 
ments (of the considered family) for which the compara- 
ble set is of this size and the question is whether there 
are some interesting differences in their performance. As 
a figure of merit for this purposes we employ the aver- 
age distance (-De), which is closely related to the quality 
of the fixed-pair comparison and partially quantifies also 
the difference of states. 

The performance of three comparison measurements 
Eswap, ^xy, and E^ is compared in Table HI It is 
clear that whichever measure /j, is used, the mean value 
{De)p is greater for the measurements Eswap and E^ 
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FIG. 4: Densities of states NsHip) (solid) and iV|f™''(p) (dashed) for different POVM effects: Sasym (left column), E^y- 
(middle column), and Ez- (right column); and for different measures: the Hilbert-Schmidt measure (top row) and the Bures 
measure (bottom row). 



TABLE I: Effectiveness of the probability-based comparison 
based on different POVM effects. Mean values (•)^ and dis- 
persions D,4-] of the distance DE{g'SiS,,S'^) for different kinds 
of measure /i. 





EswAP 




E. 


(I5e)hs 


0.07032 


0.05522 


0.07032 


{De)b 


0.09006 


0.07074 


0.09006 


Dhs [De] 


o.ofooe 


0.00695 


0.0f506 


Db [De] 


0.0f533 


0.0f062 


0.023f4 


VDhs [De]/{De)hs 


f.426 


f.510 


1.745 


VDb [De]/{De)b 


f.375 


f.457 


1.689 



than that for the measurement ^xy Furthermore, al- 
though both POVMs EswAP and lead to the same 
expectation values (-De);^, the former one gives rise to 
less dispersion [De] and relative standard deviation 
a/D^ [£)e]/(£'e)m. In addition, from Fig.|4]it follows that 
the effect iJasym results in the smallest density of states in 
the vicinity of the point po- Such feature is very demand- 
ing because the values close to pq are the most affected 
by potential statistical errors, which are unavoidable in 
practise. 

Using the mentioned figures of merit, we can draw a 
conclusion that the measurement Eswap performs (on 
average) the best among the considered measurements. 
There is yet another fact in favor for this. Fig. |4] contains 
also the densities of the same states A^|f {p) defined by 

Ap^O Ap JtT[Erj(g)ri]e[p;p+Ap] 



The value of the quantity A''|?™°(p)Ap tells us the num- 
ber of states of the form t] ^i] for which the probability 
tT[Er] <Si rj] belongs to the region [p;p + Ap] (explicit for- 
mulas for involved effects E are given in the Appendix). 
One can clearly see from Fig.|4]that for the SWAP-based 
comparison (unlike the other two) the distribution is con- 
centrated far from the border point pq. It is evident 
that if the measured experimentally probability p satis- 
fies p € P^ym then one cannot judge on states being the 
same or different. However, even in this case it is possible 
to extract an additional information. In fact, the mea- 
sured probability p sets a limit on the maximum trace 
distance between the states g and ^ (providing they are 
different), because trig — ^| < 2y^. Consequently, the 
smaller the measured probability p the closer the states 
g and ^ are. 



C. "Non-diagonal" qubit measurements 

We have argued that "diagonal" qubit measurements 
are not able to decide on the difference of states defined 
by codirectional Bloch vectors, in particular, the states 
■^I^g are not in the comparable sets of any measurement 
from this family. We addressed the question whether this 
feature is general. In other words, whether there exists a 
two-valued qubit measurement allowing us to decide on 
the difference between the complete mixture and some 
other state. 

Consider a qubit state ^ and its spectral decomposi- 
tion ^ = X]j=i,2 "^ith eigenvalues Si > S2. 
Suppose a two-valued measurement E with the effects 
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FIG. 5: States g inside the Bloch ball (determined by vectors r) which can be distinguished from a particular fixed state 
^ — diag(Hi,H2) (given by vector k = (0, 0, Hi — H2)) by using the non-diagonal POVM effect E = \E2 Hi)(H2 ® Hi| = 
diag(0,0, 1,0) for the following cases: ^ = diag(l,0) (a), ^ — diag(i, (b), ^ — diag(|, |) (c), and f = diag(j, |) (d). 



El — |ii2 ® Si)(^2 (8)^11 and E2 — I — Ei. Then we have 

De{q®£„S"^) = 2 mf |e22Si-77ii?722| 

r 2^?22Sl - i if Q22 > (4Si)-i, 

1^ otherwise, ^ ' 

where Qu = {E^\g\Ei) and ry^ = (Si|ry|Sj), i = 1,2. If, for 
instance, 2i = 1 (i.e. ^ is the north pole of the Bloch 
ball), then states g with P22 < 1/4 form the compara- 
ble set for ^. They lie below latitude 60° North and 
include also the maximally mixed state (see Fig. [5^) ■ If 
^ is the maximally mixed state (Si = S2 ~ then it 
is unambiguously distinguished from any state from the 
southern hemisphere of the Bloch ball (Fig. [SId). Ap- 
plying unitary transformations of the form UgW and 
{U (E) U)Ei^2{U^ ^ W), we can draw a conclusion that 
maximally mixed state can be effectively compared with 
any other qubit state, which answers our question in a 
positive way. 

We find that on average the fraction of comparable 
states is smaller than i. In particular, |5(Jo,„p|hs — 
1(6 - 71n2) = 0.43 and |5-^p|b = 0.42, which means 
that the measurements EgwAP, ^xy, and outperform 
the non-diagonal measurement in this parameter. The 
average distance (-De) reads values 0.1342 and 0.1524, 
respectively. We can see that the quality factor {De) is 
increased by the expense of the lower universality factor 

l^comp I ■ 

Surprisingly, for a fixed pure state ^ — |S)(S| it is 
possible to design a two-outcome measurement Eh such 
that D-B.{g (S> £.,S~^) > for all states g ^ In fact, 
suppose a POVM Eh with effects -Bhi = diag(i, |, i, |) 
and i?H2 = I — E^i specified in the orthonormal basis 
{|S ® S), |S ® |S_L ® S), |S_L «) Then, psiig® 

= 1(2-^2) andpHi(?7«'?7) = ^^ + ^"^22)^ where £122 = 
{E±\g\E±) and 7/22 = (Sj^|?y|S_L). Hence P^^ = [i, j) and 
Pij^ [ji §]■ Therefore, li g ^ S,, then we necessarily 
observe a probability pni outside the interval P~i, which 
unambiguously identifies the difference of g and ^. Such 



a two-outcome experiment can be used to check whether 
a copy g of the etalon pure state ^ was really produced. 
Nonetheless, in spite of the seeming effectiveness of this 
measurement, its average performance is quite low. To be 
precise, the fraction of comparable states g^S, on average 
reads |5(;o,np|HS = 0.097, or |5(;;„jp|B = 0.131, and the 
quality factor is {Ds)hs = 0.0049, or {De)b = 0.0079. 



V. ALMOST UNIVERSAL COMPARISON 
MEASUREMENT 

We have shown in Sec. Illll that any universal compar- 
ison measurement is necessarily an LIC measurement. 
However, the question of an existence of an almost uni- 
versal comparison measurement (which is not LIC) re- 
mains open. The maximal fraction of pairs of qubit states 
g, ^ which can be distinguished unambiguously by means 
of a two-outcome POVM seems to be i. In what follows 
we will give an example of a three- valued POVM, for 
which this fraction equals 1. This makes such a three- 
outcome measurement almost universal. 

Consider a measurement Ea.out composed of effects 
El = diag(0, 1,0,0), E2 = diag(0, 0,1,0), and E3 = 
diag(l, 0, 0, 1) in some Hilbert space basis \j ig) k) of 
two qubits. Assuming the states are the same, the 
probabilities of outcomes satisfy — < I/4, 

psamc ^ 2 - j9f'"'= - pf-""" > 1/2. In othcr words, the set 
of twin-identical states is mapped onto a line inside 
the probability simplex (see Fig. IH]). On the other hand, 
the elements of S~ , g^£,, give rise to probability vectors 
^•diff intersecting the line of twin-identical states, Pj^out: 
if and only if density matrices g and ^ have the same di- 
agonal elements. However, the subset of the states g^ £, 
satisfying this peculiar requirement has zero measure in 
. Hence, almost all pairs of states g and ^ can be 
compared by the described three-outcome measurement 
Es-out (Fig.®. 

Let us note that the considered example of 3-outcome 
POVM is nothing else but a coarse-graining of a local 
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FIG. 6: Almost universal comparison measurement (for 
qubits) with three outcomes. States q and ^ can be distin- 
guished whenever their Bloch vectors satisfy Tz 7^ fcz. 



projective measurement applied on each of the system 
independently. In particular, Ei = Eqi, E2 = Eiq, and 
E3 = Eoo + Ell, where Ejk = \j ^ k){j k\ (j, /c = 0, 1) 
are the effects forming the local (factorized) projective 
measurement E4_out- In other words, if one performs the 
same (along the same direction) Stern-Gerlach experi- 
ment on both spin-i systems, then the resulting four- 
outcome POVM E4_out performs an almost universal com- 
parison. So does the POVM Ea.out, where the outcomes 
'00' and '11' are unified into a single one. 

For a given state ^ the set of comparable states g equals 
the whole Bloch ball except for a circle of states satisfy- 
ing rz = kz, where = tr[gi(T2] and fcz — tr[<^o'z] are 
z-components of the corresponding Bloch vectors r and 
k, respectively. Here, we assume that cr^ is the local 
measurement performed on each of the states. A direct 
calculation of the distance ([T]) yields 



A-i-out((?®e,5-' 



If \rz~kz\ ifr.fc, >0, 
2 \ \rz ~ kz\ + \rzkz\ otherwise. 

(33) 

It is not hard to see that the calculation of |5j.Q^p| by 
formula ^ results in 1, i.e. the comparison measure- 
ment E3_out is almost universal indeed. The average dis- 
tance dSI) is (£'3-out)HS = 0.29, or (£>3-out)B = 0.33, i.e. 
substantially greater than for 2-valued measurements (cf. 
SecHVl). 

This qubit example can be straightforwardly general- 
ized to any dimension. Suppose a d-dimensional sys- 
tem and let Q be a projective measurement associated 
with effects Qj — where form an or- 

thonormal basis in H. Suppose the same measurement 
is performed on both d-dimensional systems and define a 



coarse-grained POVM 
Qj and Ejk = Qj ® C 
same, then 



Ecg with effects Eq = X]j=i Qj ® 
}k if j 7^ k. If the states are the 



Clearly, p 



diff 



^samc 
Pkj 



and p^^™^ = 1-2J2 



^same 
k fjk ■ 



p+ if ii^Mi^j) 



eg 



= i^^m) for all 



J = 1, . . . , d. Thus, P(+ is a less-parametric subset of 
and, hence, its measure is zero. In conclusion, the mea- 
surement Ecg consisting of 1 -I- d{d — 1) effects {Eq, Ejk} 
is an example of the almost universal comparison mea- 
surement for c?-dimensional quantum systems. 



VI. 



SUMMARY 



In its essence the comparison is a binary decision prob- 
lem, which we believe plays a very important role in 
our everyday lives. In this paper we addressed its quan- 
tum version, namely, a comparison of a pair of unknown 
sources of generally mixed states. We designed a new 
comparison strategy based on the observed statistics (not 
individual outcomes) of a particular comparison measure- 
ment device. It turns out that, basing upon the observed 
probabilities of measurement outcomes, one can some- 
times draw an unambiguous conclusion on the difference 
between the states. In fact, it seems that a vast major- 
ity of measurements are capable to compare some pairs 
of states. However, as we have shown in this paper, the 
universal comparison of two arbitrary states requires lo- 
cally informationally complete measurements, hence, a 
complete tomography of both sources is necessary and 
sufficient in order to perfectly distinguish between twin- 
identical states rj^rj and non-identical states g^S, {g ^ ^). 

Further we analyzed the comparison performance of 
two-valued qubit measurements. We defined the fam- 
ily of "diagonal" measurements including the so-called 
SWAP-based comparison measurement, which is known 
to be useful for the single-shot unambiguous compari- 
son of pure states. We have shown that none of these 
diagonal measurements is able to decide on the differ- 
ence of a completely mixed state from any other mixed 
state. Consequently, the fraction of comparable states, 
I'^complEj is at most ^ for diagonal measurements. We 
compared in detail the average performance of three diag- 
onal measurements Eswap, ^xy, and E^ that are bound- 
ary (extremal in case of Eswap and Ej) points of diagonal 
measurements. Although for all of these measurements 
I'^complE = \^ we found differences in the distribution 
of distances D^{g ® S,iS^)- Using these considerations, 
we concluded that the SWAP-based comparison performs 
(on average) better than the other two examples. 

We also provided non-diagonal comparison measure- 
ments enabling us to decide on the difference between an 
arbitrary state ^ and the complete mixture i/ or between 
a pure state ^ and any other state (in both cases the mea- 
surement depends on £). In this sense, for a given ^ the 
measurements of this kind overcome the performance of 
any diagonal measurement. However, their average per- 



formance over the set of all states result in \S^, 



E< 



-'comp I >= ^ 2 ■ 

Despite this shortcoming, any pair of qubit states can 
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be compared in a suitable non-diagonal two- valued mea- 
surement. 

In the remaining part we gave an example of the almost 
universal comparison measurement in any dimension, i.e. 
a measurement for which the size of the comparable set 
is maximal, l^^^j^pjE = 1, but still there are some pairs 
of states (forming a subset of measure zero), for which 
their difference cannot be certified. This measurement 
has d{d — 1) -I- 1 outcomes (3 in case of qubits), hence, 
it is not locally informationally complete, which requires 
the specification of 2(rf^ — 1) independent parameters (6 in 
case of qubits) and, consequently, the measurement with 
at least 2(P — 1 outcomes (7 in case of qubits). More- 
over, this measurement is experimentally very feasible to 
implement, because it is just a coarse-graining of a local 
measurement where both systems are measured by the 
same (d- valued) projective measurement (e.g., the Stern- 
Gerlach apparatus oriented along z-direction in the case 
of spin particles). 

In summary, we have shown that the universal com- 
parison of states (in the considered settings) is not pos- 
sible, but still there exist simple almost universal com- 
parators. In particular, for qubits the two-outcome mea- 
surements are not sufficient for almost universality, but 
three-outcome ones are already suflScient. A nice fea- 
ture of the proposed almost universal comparison is its 
experimental simplicity. We left many open questions, 
especially concerning an optimality of the almost univer- 
sal comparison measurements. As concerns the universal 
comparison, an optimality is directly related to the opti- 
mal complete tomography. 
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Appendix A: Densities of states 



The densities of states A^]^'^(p) and A''|f™°(p) are intro- 
duced in Eqs. and (ISTI) . respectively. It is worth not- 
ing that the domain of functions A^^'^(p) and N^^^{p) is 
Pj^ and Pj^, respectively. In this section, we present the 



explicit formulas of these densities for the effects -Easym, 
Exy-, and E^- specified in Eqs. We calculate 

the densities by using either the Hilbert-Schmidt mea- 
sure (I^B)) or the Bures measure (|27p and the obtained 
densities are depicted in Fig. |4l 

As far as POVM effect -Easym is concerned, -Pasym — 

J and the 



(0,i] and P+y m [0,i], consequently po 



density of states N!^sym,iiip) symmetrical with respect 
to the point pq. The density of states can be calculated 
explicitly in the corresponding domains for the Hilbert- 
Schmidt measure and expressed in quadratures for the 
Bures measure, namely. 



n: 



diff 

asym,HS 



n: 



diff 



asym.D 

(P) 



9 
2 

32 



1 + (4p- 1)2(2 In |4p-l|-l) 
1 



|4p-l| 



- (4p - ly 



1 



dr , 



Tvrsame 
-'*asym,HS 



6v/l-4p, 



N: 



asym,B 



(p) = A^l-Ap/ir^. 



Similarly, for the effect E^y- we have P^y^ — (0, 5], 



Pi_ = [0, i], and po 



yxy~ iiuvv. -L ^y_ — y^, 2J 

xy- — L"i4Jj m — \ but the densities of states 
differ from those obtained above and in the corresponding 
domains they read: 



^dif? 



xy-MS\ 



{p) ^ - (1 + 2(4p - l)2)arccos|4p- 1| 



-3|4p-l|Vl-(4p-l)2 



TT— arcsin |4p— 1| 



de 



arcsin |4p— 1| 
1 ' 



/r2 sin^e- {4p~ 1) 
i^p-ii y (1 - r2) sin^ 9 



dr, 



Ki'^:Bip) - 4. 

Finally, the effect Ez- is characterized by regions 
P~_ = (0,1], P^_ = [0,^], Po — ^ and gives rise to 
the following density of states: 



n: 



diff 
z-,HS 



ip) 



9 r 

4 

NT-^isip) 
NT-n{P) 



(l + (2p-l)2)(l-ln|2p-l|)-2 



V(r2-(2p-l)2)(l-r2) 



dr. 



\2p-l\ 



3p/^l-2p, 

- 2p. 
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